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Abstract 
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and the S'?7(2)-invariant ladder models are investigated. It is shown 
that corresponding to these v4 n -symmetric chain models and SU(2)- 
invariant ladder models there are exactly solvable stationary discrete- 
time (resp. continuous-time) Markov chains with transition matrices 
(resp. intensity matrices) having spectra which coincide with the ones 
of the corresponding integrable models. 
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1 Introduction 



Integrable lattice models such as quantum chain and ladder models have played significant 
roles in statistical and condensed matter physics. Many of these models can be exactly 
solved in terms of an algebraic or coordinate "Bethe Ansatz method" |IJ, see e.g., for 
chain models with periodic boundary conditions and fixed boundary conditions, and || |J 
for 2-leg ladder models with open boundary conditions. The intrinsic symmetry of these 
integrable chain models plays an essential role in finding complete sets of eigenstates of 
the systems. 

Stochastic models like stochastic reaction-diffusion models describing coagulation- 
decoagulation, birth-death processes, pair-creation and pair-annihilation of molecules on 
a chain, have attracted considerable interest due to their importance in many physical, 
chemical and biological processes f|. E.g., the simplest models for diffusion-reaction pro- 
cesses describe particles stochastically hopping on a lattice @. These diffusion- reaction 
models have been studied in various ways J7|. In particular, they have been connected to 
spin- 1/2 Heisenberg quantum spin chains ||, and then further developed and generalized 
to SU(2) symmetric spin-s chains ||. It was shown in |1(]] that the U q SU(p/m) invariant 
models Jll[ also naturally appear as time-evolution operators of chemical systems and 
the U q SU(3/0), U q SU{l/2) and U q SU(2/l) symmetric chains were discussed from a sim- 
ilar point of view in [12|]. In |13[ an exact solution of a reaction-diffusion process with 
three-site interactions (with a special next to nearest neighbour interaction) is presented. 

The stochastic reaction-diffusion systems are studied in terms of the "master equa- 
tion" which describes the time evolution of the probability distribution function JTDL IT3 



This equation has the form of a heat equation with potential (a Schrodinger equation with 
"imaginary time"). For some reaction-diffusion processes the "Hamiltonians" in the "mas- 



ter equation" coincide with the generators of the Hecke algebra []T5| . If an integrable sys- 
tem with open boundary condition can be transformed into a stochastic reaction-diffusion 
system, e.g., by a unitary transformation between their respective Hamiltonians, looked 
upon as self-adjoint operators acting in the respective Hilbert spaces, then the stochastic 
model so obtained is exactly solvable with the same energy spectrum as the one of the 
integrable system [[10], |l|, |Tj| . 

We have discussed the integrable chain models in [17] and ladder models in j3]. In this 
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paper, we give a systematic description of these models and extend the results to general 
square lattice models having a certain Lie algebra or quantum Lie algebra symmetry and 
their corresponding stationary discrete-time and continuous-time integrable stochastic 
lattice models. 

We consider M x L square lattices and give the construction of Lie algebraic (resp. 
quantum Lie algebraic) invariant lattice models in section 2 (in section 3). In section 4 we 
discuss the A n symmetric integrable chain models and the SU(2) symmetric integrable 
ladder models in the fundamental representation. In section 5 we prove that these A n 
(resp. SU(2)) symmetric integrable chain (resp. ladder) models can be transformed into 
both continuous-time and discrete-time Markov chains. Some conclusions and remarks 
are given in section 6. 



2 Lattice Systems with Lie- Algebraic Symmetry 

Let A be a bi-algebra with linear operators multiplication m and coproduct A such that 
m : A® A — > A, A : A — > A® A. Let id denote the identity transformation, id : A — > A, 
p the transposition operator, p : A <S> A — > A ® A, p(a ® b) = (b <g> a), Wa,b £ A. 
The multiplication m is associative, m(m <g> id) = m(id ® m), but not commutative in 
general, mop/m. The coproduct operator A is an algebraic homomorphism, A(ab) = 
A(a)A(6), Va, b G A. A(a) and A (6) belong to A <S> A. The multiplication of tensors is 
defined by (oi ® a 2 ) (&i ® ^2) = a ibi <8> 0.2^2, a i, a 2, 61, &2 G A The coproduct is associative 
(A ® id) A = (id (g> A) A, but in general not co- commutative poA ^ A. The operation A 
preserves all the algebraic relations of the algebra A. It gives a way to find representations 
of the algebra A in the direct product of spaces. If a bi-algebra has in addition unit, counit 
and antipode operators, it is called a Hopf algebra. Lie algebras are Hopf algebras with 
A co-commutative. Quantum algebras are Hopf algebras that are not co- commutative, 



see e.g. |T8[ and references therein. 

A Lie-algebra A is a bi-algebra. Let e = {e a }, a = 1,2, ...,n, be the basis of A, 
satisfying the Lie commutation relations 

K, e p ] = C2/3 e 7> (1) 
where C^p are the structure constants with respect to the base e. 



Let A (resp. C(e)) be the coproduct operator (resp. Casimir operator) of the algebra 
A. We have 

[C(e),e Q ] = 0, a = l,2,...,ra. (2) 
The coproduct operator action on the Lie algebra elements is given by 

Ae a = e a <g 1 + 1 <g e a , (3) 

1 stands for the identity operator. It is easy to check that 

[Ae a ,Ae f) ] = CZ fl A( h . 

From the properties of the coproduct, AC(e) is a two-fold tensor satisfying 

[AC(e),Ae a ] =0, a = 1, 2, n. (4) 

We consider M x L square lattice systems. To each point at the i-th rung, % = 1, L, 
and 0-th leg, 9 = 1, M, of the lattice we associate a (finite dimensional complex) Hilbert 
space Hf. Let — H} <g Hf <g ... <g Hf 4 . We can then associate to the whole lattice the 
tensor product Hi <g Hi <g ... ® Hi- The generators of the algebra A acting on this Hilbert 
space associated with the above lattice are given by E a = A ML ~ 1 e a , a = 1, 2, n, where 
we have defined 

A m = (1 <g> ... <g> 1 <g>A)...(l ® 1 (g> A)(l ® A)A, VmGJV. (5) 

m times 

E a also generates the Lie algebra A: [E ai Ep] = C] s pE 1 . 

Let A™ be an m-fold tensor operator with operator A on the n-th, 1 < n < m, 
tensor space and identity on the rest. For instance, A\ — A, A\ — A <g) 1, A\ — 1 <g) A, 
A\ = A <g> 1 ® 1, A^ = 1 ® A <g> 1, A^ = 1 <g> 1 ® A. Set 

M 2 1 

^ = E - E E ^...^A^-A? 2 AlC(e), (6) 

«M = 1 12 = 1 «1=1 

where a ili2 ,„ iM E € such that ft, is hermitian. Let IF denote a real entire function defined 
on the ML-fold tensor space A <g> A <g> ... (g) A of the algebra A. We call 

H=J2lF(h) t , l+ i (7) 
i=i 
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the (quantum mechanics) Hamiltonian associated with the lattice. Here IF(h)i t i+i means 
that the 2M-fold tensor element F(h) is associated with the i and 2 + 1-th rungs (columns) 
of the lattice and acts on the space Hi <g> H i+ i, i.e., 



F(h) 



Li+l — J-l 



Li-1 



F(h) ® 1 



i+2 



(8) 



where lj = 1* <g) ... ® lf^ is the identity operator acting on the i-th rungs of the lattice. 

[Theorem 1]. The Hamiltonian if is a self-adjoint operator acting in Hi (&H2 ®...®Hl 
and is invariant under the algebra A. 

[Proof]. That H is self-adjoint is immediate from the definition. To prove the invari- 
ance of H it suffices to prove [H, E a ] — 0, a = 1, 2, n. 

From the formula for the coproduct we have 



L-1 
i=l 



where (e a )j = li <g> ... ® <g> A M - 2 e a ® 
M — 1, A -1 is understood as identity operator. 

From QD it is direct to prove that [/i, A 2M ~ 2 e a ] = 0. Obviously []F(h) i)i+1 , (e a )j\ = 0, 
Vj ^ i,i + 1. Therefore we have, for all a = 1, 2, n: 



(9) 



1 L , A M ~ 2 as defined in (ED. For 



[H,E a ] 



L-1 



L-1 



^a)fc 4" ifia)i "H (Ca)j+1 



i=l 
L-1 



L-1 



(10) 



£ [^(/Om+i, (e«)< + (e a )i+i] = E [^Okm, (A 2M " 2 e 



a)i,i+l 



8=1 



8=1 



3 Lattice Models with Quantum Lie Algebraic Sym- 
metry 

Let e = {e a , h a }, a = 1, 2, n, be the Chevalley basis of a Lie algebra A with rank 
n. Let e' = {e r a , f' a , h' a }, a = l,2,...,n, be the corresponding elements of the quantum 
(q-deformed) Lie algebra A q . We denote by r a the simple roots of the Lie algebra A. The 
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quantum algebra generated by {e' a , f' a , h' a } is defined by the following relations JTB|: 

[ti a ,tip\ = 0, [h' a , e' p ] = a aP e' pi 

A h' —d h' 



q Ua — q~ 



together with the quantum Serre relations 

l-0„ 



E (-i) 7 

7=0 

1— a a g 

E (- 1 ) 7 

7=0 



1 — a a/3 
7 

1 — a a p 

7 



q a a 



where for m > n G W, 



m 
n 



m 



[n] g ![m - n] q V 



n „! = n 



["-1],-[2U1] S 



n „— n 



_ q--q 

i ~ -1 ' 

g-g 



(a Q/ g) is the Cartan matrix, 



q is a complex quantum parameter such that q da ^ ±1, 0. 

The coproduct operator A' of the quantum algebra A q is given by 



(11) 



(12) 



A% 
A'e' r 
A7, 



a 



h'®l + l®ti n 



q -d a h' a + g d a h' a e t 



f' a ® g" dah " + g dQ?i « ® /£. 



(13) 
(14) 
(15) 



It is straightforward to check that A' preserves all the algebraic relations in (|TT|) and (0). 
Let C q (e') be the Casimir operator of A q , i.e., [C g (e'), a] =0, Va G A g . For any entire 



function IF of C q (e'), we have 



[F(C ff (e')),a] = 0, Vaei, 



(16) 



and 



[A'F(C q (e% A' a] = 0, Va G A q . 



(17) 
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Especially, by formula (|I3D one gets 



£t ±d a h' a _ Q ±d a h' a gj Q ±d a h' a 



Hence 



[A'lF(C q (e')), A'q ±d ^} = [A'F(C q (e')),q 



Q 



q 



±d a h' a g q±d a h' a ^ _ q 



18) 



(19) 



[Theorem 2]. The lattice model defined by the following Hamiltonian acting in Hi 
H2 ® ... <8> -f^L is invariant under the quantum algebra 

L-l 



(20) 



2 1 



where 

^ = E - E E a ili2 .., M A'J...«C ? (e'), 

«M=1 12=1 11=1 

with aj li2 ...j M G C such that /i g is hermitian. 

[Proof]. The generators of A q on the lattice are given by 



k = A ,ML - 2 h' a = x; ii ® ... ® ® (A' 2M - 2 ^) ili+ i ® i i+2 ® ... ® ii, 



L-l 



E' 



-- q dah ' a ® ••• ® <? da/l " ® (a 

i=l 
L-l 

K = E^ Q ^®-®^®( A ' 2M " 2 ^W® ( ? 



(21) 



-d a h' a g g q—d a h' 



i=l 



From dTJ and (0) we have [fc„ A ,2M "X] = [h q , A !2M ~ 2 e' a ] = [h q , A' 2M ~ 2 f a ] = 0. 
Therefore 



L-l ji-2 L-l 

E^W. E + E ] (? 

i=l \j=l j=i+2 



d a h' a 



(A 



I2M-2J \ 

e a)j,j+l 



<g> q- dah '^ 



+q 

L-l 



d a h' 



{A' 2M ~ 2 e' a ) iti+l ®...®q- d ^ 



^2[lF(h q ),A' 2M - 2 (e> a : 



i=l 



i,i+l 



0. 



[H q ,F^\ = is obtained similarly. = can be proved like fllPp. Hence H q 

commutes with the generators of Aq. ■ 

The Hamiltonian system ( |2"(ip is expressed by the quantum algebraic generators e' = 
(h' a , e' a , /4). Assume now that e — > e'(e) is an algebraic map from A to A q (we remark that 



7 



for algebras with three generators like Ax, both classical and quantum algebraic maps can 
be discussed in terms of the two dimensional manifolds related to the algebras, see [19| ) . 
We then have 

L-X 

H q = Yl F{h g {e'{e)) i>i+ x. (22) 

In this way we obtain lattice models having quantum algebraic symmetry but expressed 
in terms of the usual Lie algebraic generators {e a } with manifest physical meanings. 



4 Integrable Lattice Models with Lie algebraic Sym- 
metry 

4.1 Quantum Yang-Baxter Equation 



The quantum Yang-Baxter equation (QYBE) |20| is the "master equation" for integrable 
models in statistical mechanics. It plays an important role in a variety of problems in 
theoretical physics such as the study of exactly solvable models like the six and eight 
vertex models in statistical mechanics [pTJ], of integrable model field theories [22j|, of exact 



S-matrix theoretical models [23~|, as well as in the investigation of two dimensional field 
theories involving fields with intermediate statistics EM, in conformal field theory and 
in the study of quantum groups [IS]. In the following we investigate the integrability of 
lattice models having a Lie algebraic symmetry constructed in section 2. We also present 
a series of solutions of the QYBE from the construction of integrable models. 

Let V be a complex vector space and R the solution of QYBE without spectral pa- 



rameters, see e.g. [[18]. Then R takes values in End@(V ® V). The QYBE is 



-R12-R13-R23 



-^23-^13-^12- 



(23) 



Here Rij denotes the matrix on the complex vector space V <8> V ® V, acting as R on the 
z-th and the j-th components and as the identity on the other components. 

Let R = Rp (p is the transposition operator). Then the QYBE (|23|) becomes 



-R12-R23-R12 



R23R12R23, 



(24) 



where Rx2 — R ® ly, R 



23 



R and ly is the identity operator on V. 
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In the following we say that a lattice model with nearest neighbours interactions having 
a (quantum mechanical) Hamiltonian of the form 

H = sWi+i (25) 

i=l 



is integrable in the sense that the operator 7i satisfies the QYBE relation (EM), i.e 



(ft)i2(ft) 23 (W)i2 = {H) 2 3{H) 12 {H) 23 , (26) 

where (TC)i2 = Ti- <8> ly an d (^)23 = ly ® Here 7i is a solution of the Yang-Baxter 
equation without spectral parameters. After "baxterization" the Hamiltonian system 



(|25|) satisfying relation (|26|) can in principle be exactly solved in terms of algebraic Bethe 
Ansatz method, see e.g. Here the vector space V is taken to be the Hilbert spaces 
associated with one rung of the lattice. 



4.2 Integrable A n Symmetric Chain Models 

The integrability of the models having a Lie algebraic symmetry presented in section 2 
depends on the detailed representation of the corresponding symmetry algebra. In this 
section we investigate the integrability of chain models with nearest neighbours interac- 
tions and Lie algebraic symmetry A n . 

Let (a a p) be the Cartan matrix of the A n algebra. In the Chevalley basis the algebra 
A n is spanned by the generators {h a , e a , f a }, a = 1,2, ...,n, with the following algebraic 
relations: 

[h a , hp] = 0, [h a , ep\ = a a p€p, [h a , fp\ = —a a pfp, [e a , fp\ = S a/ 3h a , (27) 
together with the generators with respect to non simple roots, 

e Q .../3 7 = [e a , [ep, e 7 ]...], fa.../3-y = [fa, [f(3, / 7 ]---]- (28) 

Let E a p be an {n + 1) X (n + 1) matrix such that (E a p) r s = ba-fbps, i-e., the only non 
zero element of the matrix E a p is 1 at row a and column j3. Hence 

E Q pE 7 s = Sp^E a s, [E a p, E^s] = 5p~/E a s — 5s a Ep 7 . (29) 
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For the fundamental representation we take the basis of the algebra A n as 

h a = E aa — E a+ i ia+ \, a — 1, 2, n 

e = {E a(3 } ) (30) 

> (3 > a— 1,2, ...,n 
f = {E Pa } J 

Both {e a } and {f a } have a total of n{n + l)/2 generators. 

With respect to the basis (|30|), the Casimir operator of the algebra A n is given by 

n(n+l)/2 n 



a=l a=l 



C An = (n + 1) 51 ( e «/« + Z" e «) + H a ( n + 1 _ a )^o 



n n— a 



(31) 



+ X! X! 2a(n + 1 — a — f3)h a h a+ p — a, 
a=l p=\ 

where a is an arbitrary real constant. 
The coproduct operator A is given by 

A(1) = 1®1, %) = /i a ®l + l®/i a , a =1,2,..., n 

A(e /3 ) = ep®l + l®ep ) (32) 

/3 = l,2,..,n(n + l)/2, 

A(/^) = ^(8)1 + 1®/^ J 
where the identity operator 1 is the (n + 1) x (n + 1) identity matrix. 
By (|31|) and © we have 
AC An = C An ® 1 + 1 ® C A „ - al <g> 1 



n(n+l)/2 „ 

+(n + 1) 51 ( e « ® /« + /« ® e Q ) + a ( n + 1 - a )^a ® 

?i n— a 

+ X! H "( n + 1 - a - /3)(^a ® + ^a+/9 ® K) ■ 

a=l /3=1 



(33) 



It is easy to check that under the representation ( |30D C An is equal to n{n + 2)1. 
Therefore the sum of the first two terms on the right hand side of (|33|) is 2n(n + 2)1 x 1. 
In the following we take a in (^) to be 2n(n + 2) so that the terms that are proportional 
to the (n + l) 2 x (n + l) 2 identity matrix will disappear in (|33|). 
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From ( |5DD and ( |3"3"| ) we have 

n+l 

AC An = (n + l) ]T £ a/3 <g> ^ 
n 

+ ^ a(n + 1 - a)(E aa - E a+1<a+1 ) ® (£ aQ - £ a+ i )Q:+1 ) 

n n— a 

+ ^2^2 a ( U + 1 — « — /5)[(-E QQ — E a+ i tCC+1 ) <g> (£' Q , +J g iQ , +( g — E a+/ 3 + i )Q , +/ g + i) 
a=l /3=1 

+ (-E'o+/3,a+/3 — ^a+/3+l,Q!+/3+l) ® (-^aa — -^a+l,a+l)]- 

(34) 

AC J 4„ in (0) is an (n + l) 2 x [n + l) 2 matrix. Its matrix representation is 

(AC^Ja/3 = S a p[(n + l)5 a ,l(n+l)+l+l ~ 1] 

+ (n+ l)[8 at j^ n+ 2)+k+2^,(J+l)(n+2)+k(n+l) (35) 
+^,j(n+2)+fc+2^a,(j+l)(n+2)+fe(n+l)], 

where a, [3 = 1, 2, {n + l) 2 , / = 0, 1, n, j = 0, 1, n — 1, k = 0, 1, n — j — 1, 
8 a ,j{n+2)+k+2 — if a 7^ j(n + 2) + k + 2 for all possible values of j and k. For example, 

/ 1 \ 
0-120 
2-10 
V 1 / 

[Lemma 1]. AC An satisfies the following relation 

(AC An ) 2 + 2AC An -n{n + 2)l®l = 0. (37) 

[Proof]. From (F3^) we have 

(n+l) 2 

[(AC A J 2 ] Q7 = £ (AOO^ACaJ^ 

,3=1 

= 5 Q7 [(?2+1) 5 a ,«(n+l)+i+l57,i'(n+l)+i'+l] 

— (n + 1) (5 a ,l(n+l)+l+l + &y,/'(n+l)+2'+l) + 1] 

-2(n + l)[5 

a ■ 

+(n + l) 2 [5 a ,j(n+2)+A:+257,j(n.+2)+fc+2 a. (i+l)(n+2)+fc(n+l)^7,(j+l)(n+2)+fc(n+l)] 

= -2(AC An )a 7 + in + l) 2 5a 7 [^ 

a, j(n+2)+fc+2 + S a ,(j+l)(n+2)+k{n+l)\ 



AC Al 



(36) 



+(n + l) 2 5 a7 <5 



ary 



-2(AC A J a7 + n(n + 2)5, 
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where the identity 

3a,l(n+l)+l+l ~l" 3a,j(n+2)+k+2 ~l~ ^a,(j+l)(n+2)+fc(n+l) = 1) 

(38) 

I — 0, 1, n, j — 0, 1, n — 1, k — 0, 1, n — j — 1, has been used. ■ 
[Lemma 2]. The coproduct of the A n Casimir operator AC An has the following prop- 

{AC An ®l){l®AC An ){AC An ®l) 

-n[(l ® AC A J(ACU„ ® 1) + (ACU„ ® 1)(1 ® ACUJ] (39) 
+ (n 2 - l)(AC An ® 1) + n 2 (l ® AC A J + n(l - n 2 )l ® 1 ® 1 = 

and 

(l®AC An )(AC An ®l)(l®AC An ) 

-n[(l ® AC A J(AC An <g> 1) + (AC Ar[ <g> 1)(1 ® ACUJ] (40) 
+(n 2 - 1)(1 ® AC An ) + n 2 (AC An ® 1) + n(l - n 2 )l ® 1 ® 1 = 0. 
[Proof]. By using the representation of AC An in (^) we have 

a— 7)/(n+l)+l,(/3— 7)/(n+l)+l 
= 5 aj g[(n + l)5 a _ 7) i( n+ l)( n+ 2) — 1] 

+ (n + l)[5 a _ 7i ( n+ i)(j( n+ 2)+fc+l))5/3_ 7i („ + l)(j(„ + 2)+(fc+l)(n+l)) 
+^/3-7,(n+l)0'(n+2)+fe+l)^a-7,(n+l)0'(n+2)+(fe+l)(n+l))] 

and 

(1 <g> ACA„) a /3 = (AC , A, l ) a -(„+i)2( 7 '_i) ij a_(„ + i)2( 7 /_i) 

= 5 a/3 [(n+ l)5 Q -(„ + i)2(y_i)j( n+ i) + i + i) — 1] 

+ (n + 1) [5 Q! _( Tl+ i)2( 7 /_i) ) j(„ +2 )+fc+2))^/3-(n+l) 2 ( 7 '-l),(j+l)(ri+2)+fc(n+l) 
+^-(n+l) 2 (7'-l)j(™+2)+fc+2^a-(n+l) 2 (7'-l),(i+l)(n+2)+A:(n+l)], 

(42) 

where a, /3 = 1, (n + l) 3 , I — 0, 1, n, j — 0, 1, n — 1 and k — 0, 1, n — j — 1 as 
in formula fl35|), 7 = 1, n + 1 such that (a — j)/(n+ 1) and (/3 — 7) j{n + 1) in (p|) are 
integers and 7' = 1, n + 1 in 



(41) 



Using the formulae ( |4"T]) and ( [4*2"D one can get (^) and ( |4T]| ) from straightforward 
calculations. ■ 

From Theorem 1 we know that the following Hamiltonian is invariant under A n 

H=J2lF(AC An ) iii+1 . (43) 



i=l 
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For the given representation ( |5UD of A n the integrability of Q55D depends on the form of 
the entire function IF. Due to the relation (|37D in Lemma 1, (AC^)', Z > 2, can be 
expressed as cAC^ + c'l <E> 1 for some real constants c and c' . Therefore IF(AC An ) is a 
polynomial in ACa„ up to powers of order one. 

[Theorem 3]. The following A n invariant Hamiltonian is integrable 



8=1 1=1 

L-l T n(n+l)/2 n 

(n + 1) 51 (( e a)i(/a)i+i + (/«)i(e a )<+i) + a ( n + 1 - ot)(h a )i(h a ) i+1 



E 

i=i 



a=l 



a=l 



n n— a 



+ EE a ( n + 1 - a - /3)((/i a )i(/i a+;3 ) i+ i + (/i a+ ^)j(/i a )j+i) 

a=l /3=1 



L-l, 



(44) 

where TC = AC An + 1 and the number 1 should be understood as the identity operator, 
1 <S> 1, on the tensor space Hi <g> ... ® if^. 

[Proof]. What we have to prove is that H satisfies the QYBE i.e., (H) 12 (H) 23(H) n 
(H) 23 {H) 12 {H) 23 , where (W) 12 = (AC An + 1) ® 1 and (ft) 23 = 1 ® (AC 4 + 1). We have 

(H) l2 (H) 23 (H) 12 = (1® AC A J(AC A „(8)l) + (AC^(8)l)(AC An (8)l) 

+(AC An ® 1)(1 ® AC An )(AC An ® 1) + (AC An ® 1)(1 ® AC A J 

+2AC An ® 1 + 1 ® AC An + 1 <g> 1 ® 1 

and 

(^23(^)12(^)23 = (l®AC An )(l®AC An ) + (l®AC An )(AC An ®l) 

+(1 ® ACUj(AC An ® 1)(1 ® AC A J + (AC An ® 1)(1 ® AC A J 
+2(1 ® AC An ) + (AC A „ ® 1) + 1 O 1 ® 1. 

Hence 

(W)i 2 (W)23(W)i2 - (H) 23 (H) 12 (H) 23 = 1 + 11 + III, (45) 

where 

/ = (AC A „®l)(AC A „®l)-(l®AC A J(l«)AC An ), 
J/ = (AC An ®l)(l® ACU B )(AC An ®l)-(l® ACUJ(ACU„®1)(1® AC A J, 
/// = AC An ® 1 - 1 ® AC An . 
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Using fl3"7| ) we have 

/ = (AC An ) 2 ® 1 - 1 ® (AC A J 2 

= -(2AC An - n(n + 2)1 ® 1) ® 1 + 1 ® (2AC An - n(n + 2)1 ® 1) 
= -2(ACU n ®l-l® AC A J. 

Hence 

I + 111 = l®AC An -AC An (8)1. (46) 

By Lemma 2 we get 

II = AC An ® 1 - 1 ® AC An . 

Therefore 

(H) 12 (H) 23 (H) 12 - (H) 23 (H) 12 (H) 23 = 1 + 11 + 111 = 0. 



Related to the integrable chain model ([441), there is a Temperley-Lieb (TL) algebraic 
structure in the sense that the model gives a representation of the TL algebra. An 
(L — l)-state TL algebra is described by the elements e», i = 1, 2, L — 1, satisfying the 



TL algebraic relations 25 



eie i± iei = e { , ejej = , if |i - j\ > 2 , (47) 

and 

el = fa , (48) 

where /? is a complex constant and i = 1, 2, • ■ • , L — 1. 

We suppose that the representation of an (L — l)-state TL algebra on an L chain is 
of the following form, 

tH = l x ® 1 2 ® • • • ® li-i ® -E ® lj+2 ® • • • ® U , (49) 



where 1 is the (n + 1) x (n+1) identity matrix as in section 3.2 and E is a (n + 1) 2 x 



2 



matrix. According to formulae ( [48| ) and (47) should satisfy 

E 2 =(3E. (50) 
(-B ® 1)(1 ® £)(£ ® 1) = £? ® 1, 

(51) 

(1 ®E){E® 1)(1 ® i?) = 1®£. 
14 



From the representations of the TL algebra one can construct integrable chain models 
(for the construction of the TL algebraic representations associated with the quantum 



Ai, B n , C n and D n algebras, see |26|). It is straightforward to check that for a given 
representation of the TL algebra of the form ([49]) with E satisfying (^) and (|5l]) , 

a=B+= g^gE3 181 



is a solution of the QYBE (l24f) . However in general for a given solution R of the QYBE 
fl2~4|) , there does not necessarily exist a TL algebraic representation of the form (|49|) with 
E = aR + b satisfying (50) and ([5l|) for any constants a and b. Nevertheless the solutions 
Ti of the QYBE in our A n symmetric integrable model (|44]) do give rise to TL algebraic 
representations in the following sense: 

[Theorem 4]. The following [n + l) 2 x {n + l) 2 matrix 

T-/ 

E = + 1 <g> 1 (52) 

n + 1 

gives the (L — l)-state TL algebraic representation (^) with (3 = 2. 

[Proof]. What we should check is that E in ( |52"D satisfies equations ( ^H| ) and (^Tj). By 
Lemma 1 we have 



n + 1 (n + 1) 2 



i.e., /? = 2. 

From Lemma 1 and fl39|) in Lemma 2 we get 

1) 

j ■[(AC An ®l)(l^AC An )(AC An ®l) 



(71+1) 

-«((AC An ® 1)(1 ® AC A J + (1 ® AC A J(AC An ® 1)) 
+2n(n + l)AC An ® 1 + n 2 l ® AC An - 2(n 3 + n 2 )l g> 1 g> 1] 

, ~ iVt [(w + l) 2 AC An <g> 1 - n(n + 1) 2 1 <8> 1 ® 1] = E <g> 1. 
(n + 1 r 
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By using Lemma 1 and formula ( PD| ) in Lemma 2 we conclude that 

1)(1®E) 

= T-TTV3^ ACa - ® !)(! ® ^C An )(AC An ® 1) 
(n + l)' 3 

-«((AC A „ ® 1)(1 ® AC A J + (1 ® AC An )(AC An ® 1)) 
+2n(n + 1)1 <g> AC An + n 2 AC An <g> 1 - 2(n 3 + n 2 )l g> 1 g> 1] 

= , = -rr ? [(w + 1) 2 1 ® AC An - n(n + 1) 2 1 ® 1 ® 1] = 1 ® £. 
(n + iy 

■ 

From (|52|) we see that the Hamiltonian of the A n symmetric integrable chain model 
( pH| ) can be expressed by the TL algebraic elements 

L-1 L-1 

H An = E (WW = E (" + + (n + 1)(L - 1), (53) 
j=i j=i 

with e, as in ( pE5| ) and i£ as in (0). Hence instead of the algebraic Bethe Ansatz method, 
the energy spectrum of H An can also be studied by using the properties of the TL algebra 
[p7| ] (for the case of Heisenberg spin chain model, n — 1, see 

4.3 Integrable S't/(2)-Symmetric Ladder Models 

We consider ladder models (M = 2) with SU(2) symmetry. Let Si, i = 1,2, 3, and 
C be the generators of the algebra SU(2) and the Casimir operator respectively. The 
coproduct of the algebra is given by ASj = 1 <8> Si + Si <S> 1, i = 1,2,3. Accounting to 
that A^>(e) = F{A{e), i = 1,2,3, j = 1,2, Ve e 5C/(2), the generic h is of the form 
FiC^C^Ca), where 

3 

Ci = E(s* ® 1 ® 1 ® s< + 1 ® s< ® 1 ® •S'i + 1 ® 1 ® s'i ® si), 
(=i 

3 

c 2 = E(Si ® i ® i ® Si + ^ ® si ® i ® i + Si ® i ® Sj ® i), 

i=l 
3 

C 3 = E(Si ^lOl^Si + loSi^loSi + S^loSi^l + loSi^SiOl). 

i=l 

In the spin-i representation of the algebra SU{2), the solutions of the QYBE (^) are 
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16 x 16 matrices. For instance, it is direct to check that 

n, 108d-55/ /-y i -72d+104/ n , -486(2+211/.^ , -756(2+370/ /-y 
n ° - 108 ° m H 288 ° 112 H 270 ° 113 H 216 ° 121 

-fi C 122 + ^C 123 + *=tC m + ^FC 132 (54) 

i -108(2+43/ | -216(2+80/^ , U//i , 216(2-119/ 
+ 540 ° 133 + 864 °211 + iog<-212 + Tqs ^213 

satisfies (|26|) for all d, f £ tt, where Cy*. = Cj ■ Cj ■ Ck, k — 1, 2, 3. 

The corresponding solution related to the S'?7(2)-symmetric integrable ladder model 
in || can also be expressed in the form F(Ci, C 2 , C3), i.e., 

ni 5 ^» 61/-» 41 //m /m \ , 21 n 

ft — — Ig^lll ~~ 32 < -"'112 — 30^113 ~ 48^121 — <~422j + 16 123 

(55) 

<3r< 17 n _L 173/0 _|_ 55^-y 5n , 131 

+ 4^131 — 12°132 + 240° 133 96° 211 — 3^212 + ^g-^213- 

Through baxterization, TC(x) = (x — 1)H + 16Ii6xi6 satisfies the QYBE with spectral 
parameters: Hi 2 {x)H 23 {xy)Hi 2 {y) = H 23 {y)H 12 {xy)H 23 {x) , where H 12 {-) = H{-) <g> I 4x4 , 
7^23 (•) = 14x4 ® ^( - )) I^xn denotes the n x n identity matrix. The model can be exactly 
solved using algebraic Bethe Ansatz method. It describes a periodic spin ladder system 
with both isotropic exchange interactions and biquadratic interactions: 

H = \ E(| + 2S M ■ S w )d + 2S 2 , ■ S 2 , m ) - \ £(I + 2S M ■ S 2 , +1 )d + 2S 2 , ■ S L 

i=l i=l " 

+ c E (o + ^Si j • S 2 j)(- + 2Si i+i • S 2 ,i+i), 
6^2 2 

where S^j = (ag i ,ag i ,ag i )/2, a x ,a y ,a z are Pauli matrices. Si^ (resp. S 2iJ ) is the spin 
operator on the first (resp. second) leg of the i-th rung of the ladder. 

It is also easy to see that for a more general form of 



(56) 



nil -45+23 a-4 6-28 cn i -99-3a-3b-cn , -1098-91 a-118 6-16 c n 
n - 432 Oui + 2^ On2 + ^0 °U3 

i -369-97 a-70 6+50 c ri , 396+4 a+316+25c ,^ j_ 189+29 0+20 6-4 ^ 
H 432 ° 121 H 432 0122 H 144 ° 12 3 

, 3n _l -306-2 a-29 b-\4c n i 1557-71 a+172 6+124 c n 
+ 4^131 H °132 H 2160 ° 133 

i 495-a+53 6+47c ^> , -720-22 a-49 6-43 c n i 1179+91 a+118 6+16 c n 
+ 864 ° 2 H + 432 ° 2 1 2 + 432 ° 2 13 

L-l 

with a,b,c G M, the corresponding ladder model H' = V] ■ i+1 can also be exactly 



solved by the coordinate Bethe ansatz 
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5 Integrable Models and Stationary Markov Chains 



5.1 Stationary Markov Chains 

Let us first briefly recall some concepts of the theory of Markov chains (for a detailed 
mathematical description of Markov chains, we refer to p9| ). Let Q denote the sample 
space (the set of all possible outcomes of an experiment), which we assume to be finite or 
countably infinite. Let P be a probability measure defined on the cx-algebra of all subsets 
of Q. Thus any subset A of Q is an event with probability P(A). 

Any function X = X(u), to G O, that maps the sample space into the real numbers 
is then a random variable. A stochastic process is a family (X t )tei, I a certain index set, 
of random variables defined on some sample space Q. If I is countable, i.e., I G IN, the 
process is denoted by Ai, A 2 , ... and called a discrete-time process. If / = M + , then the 
process is denoted by {X t } t > and called a continuous-time process. 

The range of X (a subset of real numbers) is called the state space. In what follows 
we consider the case where the state space S is countable or finite. In this case the related 
stochastic process is called a (stochastic or random) chain. 

Let E,F be two subsets of Q. We denote by P(E\F) the (conditional) probability of 
E given that F has occurred. A discrete-time stochastic process {A^}, % = 1,2,... with 
state space S = IN is said to satisfy the Markov property if for every I and all states 
ii,i2, H it is (a.s.) true that 

P[Xi = klX^x = ii^Xi-2 = h-2, Xi = it] = P[Xi = ii\Xi_t = h-i], 

i.e., the values of Xi_ 2 , Xt in no way affect the value of Xi, given the value of X^f 
Such a discrete-time process is called a Markov chain. It is said to be stationary if the 
probability of going from one state to another is independent of the time at which the 
transition is being made. That is, for all states i and j, 

P[X t = j\X^t =i] = P[X l+k = j\X l+k _t = i] 

for k = — (I — 1), — (I — 2), — 1, 0, 1, 2, .... In this case we set = P[Xi = j\Xi_i = i] 
and call Pij the transition probability for going from state i to j. 

For a discrete time stationary Markov chain {Aj}, i G IV, with a finite state space 
S = {1, 2, 3, m}, there are m 2 transition probabilities {pij}, i,j = 1, 2, m. P = (p^) 
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is called the transition matrix corresponding to the discrete-time stationary Markov chain 
{Xi\. The transition matrix P has the following properties: 

m 

Pij > 0, J^Pij = !> hi = 1,2, ...,m. (57) 
i=l 

Any square matrix that satisfies condition (|5"?D is called a stochastic matrix. 

A continuous-time stochastic process, {X t }t^iR + is said to satisfy the Markov property 
if for all times to < £i < ... < t\ < t and for all I it is true that 

P[X t = j\X t0 = i , X h = i u X tl = z,] = P[X t = j\X tl = k\. 

Such a process is called a continuous-time Markov chain. It is said to be stationary if for 
every i and j the transition function, P[X t+ h = j\X t = i], is independent of t. In this 
case P(t) = P(X t= j\X — i) is a semigroup (e.g. on l 2 (S)), called transition semigroup 
associated with the Markov chain. Its generator Q = (g^) has the properties: 

ftj>0, i^j, qii = -J2liji ^J = l,2,...,m (58) 

and is called an intensity matrix. Vice versa, any Q (satisfying (^) and properly defined 
as a closed operator when 5* is infinite) gives rise to a unique continuous-time transition 
semigroup, P(t) = e^', t > 0, which can be interpreted as transition semigroup associated 
to a certain Markov chain (with state space S) [^9] . 

The properties of Markov chains are determined by the transition matrix P resp. 
intensity matrix Q. If the eigenvalues and eigenstates of P resp. Q are known, then 
exact results related to the stochastic processes, such as time-dependent averages and 
correlations, can be obtained. 

Now we consider a square lattice (in the algebraic sense of sections 2-4) with ML sites. 
To every site at the z-th rung and ^-th leg of the lattice we associate N states described 
by the variable t^q taking iV integer values, 



Ti ft = T A 



-0 i 7& = 1 i 7& = 2 i ... j t3 = JV). (59) 



We associate to any lattice site a Hilbert space of dimension iV (dim H® = N). The state 
space of the algebraic lattice is then finite and has a total of (N) ML states. 

For a given integrable lattice model with Hamiltonian H, the model remains integrable 
if one adds to if a constant term c and multiplies if by a constant factor d . Moreover 
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the eigenvalues of H will not be changed if one changes the local basis, i.e., the following 
Hamiltonian, defined by 

BHB-\ (60) 

where B = ®^£fPj, B^ = B and B is an N x N non singular matrix, has the same 
eigenvalues as H. Therefore if an integrable lattice model with Hamiltonian H can be 
transformed by B (modulo constants c, d in C) into a matrix M = P resp. Q, in the 
sense that 

M = B(cH + cH)B-\ (61) 

where 1 is the N ML x (N) ML identity matrix, B as in (^), and P resp. Q as in (|57| ) resp. 
fl58|) (with m = (N) ML ), then P resp. Q defines a discrete-time resp. continuous-time 
Markov chain and the related stochastic process can be studied by using the properties 
of the corresponding integrable model with Hamiltonian H. 

In the following we discuss the question of whether the integrable lattice models ob- 
tained in the way presented in this paper can be transformed into some stationary Markov 
chains through transformations of the forms (0). We remark that in general P and Q 
have different spectra, hence it is necessary to discuss the cases M = P and M = Q 
separately. 



5.2 Discrete-time Markov Chains Related to A n Symmetric 
Integrable Models 

We first note that for an integrable chain model with Hamiltonian H = J2f=i hi^+i and 
(n + 1) states at every site i, i = 1, 2, L, if the sum of the elements in any row of the 
(n+l) 2 x(n + l) 2 matrix h is 1/(L — 1), then the sum of the elements in any row of the 
matrix H is 1. Hence if under the following transformation h — > h! given by 

ti = (B ® B){c'h + cl ® l)(B~ l <g> iT 1 ), (62) 

the sum of the elements in any row of hi is \ j[L — 1) and {h') a p > 0, a, f3 = 1, 2, (n + 
l) 2 , for some real constants c', c and a non singular (n+ 1) x (n+ 1) matrix £>, then P = 
X^^Ti 1 ^ j+i defines a stationary discrete-time Markov chain. P has the same eigenvalue 
spectrum (shifted by a constant) as the spectrum of the integrable model with Hamiltonian 
H. If P is invariant under a certain algebra A, we call the Markov chain A-symmetric. 
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[Theorem 5]. The following matrix 

(L-l)(n + l) An (L-l)(n + l)^ 



L-l 



J2 (AC An + 1 ® 1),, 



i+1 



1 



L-l 

E 



(L-l)(n + l) £j 
+ ^ a(n + 1 - Qt)(ft, a )i(/i a )i + i 



ra(n+l)/2 

n + l ) E (( e «)i(/«)i+l + (/a)t(e a )i+ly 



a=l 



a=l 
n n—OL 



a=l 0=1 

defines a stationary discrete-time A n symmetric Markov chain. 



+ 



[Proof]. I. Set ti 



i 

(L-l)(n+l) 



(AC An + 1 ® 1). Then 

L-l 

P A„ = E K,i+l' 

i=l 



From formula (EE) we have 



00* 



(L-l)(n + l) 
1 



{AC An + 1 ® 1) 



a/3 



[#0,9 [(« + l)5a,/(n+l)+Z+l - 1] 



(L-l)(n + l) 

+ (n + l)[5 a j( n+ 2)+fc+2^/3,(i+l)(n+2)+fc(n+l) 
+ ( ^/3j(«+2)+fc+2 ( 5«,(j+l)(n+2)+fc(n+l)] + ^a/S] 

+ 5aj(n+2)+fc+2^,( i+ l)(„+2)+fe(n+l) 
+^/3j'(™+2)+fc+2^a,(j+l)(n+2)+fc(ra+l)] > 0. 

Therefore (P A J a/3 > 0, a, /3 = 1, 2, .., (n + l) 2 . 
II. By using the identity (|38|), we get 



(n+l) 2 

E flu 



i 



(n+l) 2 

(L - l)(n +T) 5 [(n + ^M^iHH-a 

+ (n + l)[^aj( n +2)+fc+2^/3,(j+l)(n+2)+fc(n+l) 

+^/3,j(n+2)+A:+2^Q,(i+l)(n+2)+A:(n+l)]] 
1 1 



(L-l)(n+l) 



{n + l) 



L-l 
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Hence the sum of the elements of any row of the matrix P^ n is one, i.e., Ylp=i (PA n )af3 — 

2^/3=1 \l^i=l "'i,i+l)af3 — -L- 

III. As Ha„ is invariant under A n , Pa„ = p^nfe+i) * s obviously invariant under A n 
and has the same spectrum as Ha u - 

By the definition (|57|) P_A n is the transition matrix of a stationary discrete-time A n 
symmetric Markov chain. ■ 

The state space of this stationary discrete-time A n symmetric Markov chain associated 
with the stochastic matrix Pa„ is S — (1, 2, (n + which corresponds to {n + 1) L 
states, 

(t ® Ti ® ... <g) T n ), (66) 

Tj = ! as in (|59|) with = 1, of the algebraic chain with L lattice sites. This stationary 
discrete-time A n symmetric Markov chain model describes a chain with L sites and n + 1 
possible states, say, r/, rf + , at site i, i = 1, L. By calculation the allowed dynamics 
for any nearest neighbour pair is the interchange of their states: (r J a ,r i ^ 1 ) — > (T"f , r^-J, 
a,j3 & {l,...,n + 1}. For n = 1, there are two possible states at every site: empty or 
occupied by one particle. The A\ Markov chain describes then stochastic hopping of 
particle into left and right vacancies, which is the well known SU(2) random chain. For 
n = 2, i.e. the SU(3) case, there are three possible states at every site, say, empty, one 
spin up and one spin down particle states. This model describes stochastic hopping of 
spin up and down particles into left and right vacancies, interchanging of one spin up and 
one spin down particle at the nearest neighbours. 

The properties of a Markov chain are determined by the transition matrix P = (pij). 
A subset C of the state space S is called closed if = for all i £ C and j ^ C. If a 
closed set consists of a single state, then that state is called an absorbing state. A Markov 
chain is called irreducible if there exists no nonempty closed set other than 5* itself. A 
non irreducible Markov chain is said to be reducible. 

From formula (^) we have 

{h') aa = (/i')(n+l) 2 ,(n+l) 2 = T3T, (h') a/3 = {h')p a = 0, f3 ^ a, 

(67) 

a = l{n + 1) + I + 1, / = 0, l,...,n. 
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Let 

c (\ , (n + l)((n+ l) L ^-l)+n \ 

Do = a\a = I h 1 , I — 0, 1, n, (68) 

\ n J 

be a subset of the state space S. From formula (0), with 

(h%+i = li <g> 1 2 ® ... <g> <g> /i' ® l l+2 ® ••• ® 

we get 

(P A J aQ = l, (P A J^ a = {P An ) a p = 0, /^«, aeS . (69) 

Therefore the n + 1 states in Sq are absorbing states of the Markov chain Pa„- This 
chain is by definition reducible. For a reducible Markov chain the "long time" probability 
distribution, if it exists, may depend on the initial conditions, i.e., limj_ +00 (P J 4 n )'o may 
depend on 7. From the properties ( |69"D of Pa„, we see that if the Markov chain Pa„ 
is initially in one of the states a e So, it will remain in that state a forever. These 
n + 1 absorbing states correspond to the states of the algebraic chain through (|6"6j). For 
instance, the states 1 and (n + 1) L in S correspond to the states (0,0, ...,0) (all the sites 
of the algebraic chain are at state 0) and (n, n, ...,n) (all the sites of the algebraic chain 
are at state n). 

5.3 Continuous-time Markov Chains Related to A n Symmetric 
Integrable Models 

For an integrable chain model with Hamiltonian H = J2i=i fh,i+i an d with (n + 1) states 
at every site of the chain, if the sum of the elements in any column of the matrix h is 
0, the sum of the elements in any column of the matrix H is also 0. Hence if under the 
following transformation h — > h" with: 

h" = {B® B){c'h + cl ® lXiT 1 ® B- 1 ), (70) 

the sum of the elements in any column of h" is and (h") a ,/3 > 0, a 7^ (3 = 1, 2, (n + 1) 2 , 
for some real constants c', c and a non singular (n + 1) x (n + 1) matrix B, then Q = 
Y^i=\ h'li+i i s ^ ne intensity matrix for some stationary continuous-time Markov chain. 
Q has the same eigenvalue spectrum (shifted by a constant) as the spectrum of the 
Hamiltonian H . We call the Markov chain A symmetric if Q is invariant under the 
algebra A. 
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[Theorem 6]. The following matrix Q is the intensity matrix of a stationary continuous- 
time Markov chain: 

L-1 

Qa„ = H An -{n + l){L-l) = y £ l (AC An -nl®l) i>i + 1 

i=i 



L-1 

E 

i=l 



n(n+l)/2 „ 

(n + 1) 51 (( e «)i(/«)i+i + (fa)i(e a ) i+ i) + E a ( n + 1 ~ a ) ( h a ) i(h a ) i+1 



a=l a=l 



n n—a 



+ E E a ( n + 1 - « - (3)((h a )i(h a+ p) i+1 + (/i a+(3 )j(/i a ) i+ i) 

a=l 0=1 



n(L-l). 

(71) 



[Proof]. Set h" = (AC An - nl ® 1). Then 

Q^ = XXi+i- (72) 



i=l 

From ( p5| ) we observe that, for a ^ (3, 

= ( n + l)[<W(n+2)+fc+2<Wj+l)(n+2)+fc(n+l) 
+^j9,i(n+2)+Jfe+2^a,Cj+l)(n+2)+Jk(Ti+l)] — 0- 

Therefore (Qa„ W/3 > 0, a, /? = 1, 2, (n + 1) L . 

Again by (|35|) the sum of the elements in any given column j3 of the matrix h" is 

(n+l) 2 (n+1) 2 

E = E + l)[*o/3(<W(n+l)+i+l - 1) 

a=l a=l 

+ (n + l)[(5 aj j( n+ 2)+A;+2^/3,(j+l)(n+2)+fe(n+l) + ^/3,j(n+2)+k+2^a,(J+l)(n+2)+k(n+l)] 
(n+1) 2 

= ^ (n + l)[—8 a p + (5 aj j( n+ 2)+A;+2^,(j+l)(n+2)+fe(n+l) 

a^Z(n+l)+Z+l 

+<5/3,j(n+2)+fc+2<5a,(j+l)(n+2)+fc(n+l)] 
= (n + 1)( — 1 + <5/3,(j+l)(n+2)+fc(n+l)U=j(n+2)+fc+2 + <5/3,j(n+2)+fc+2 1 a=(j+l)(n+2)+fc(n+l) ) 

= o, 

i.e., the sum of the elements in any given column (3, (3 = 1, 2, (n + l) 2 , of the matrix h" 
is zero. Therefore the sum of the elements in any given column (3, (3 = 1, 2, (n + 1) L , 
of the matrix Q An is also zero, E^Ji^ (QaJoP = 0. At last Q An = H An - (n + 1)(L - 1) 
is obviously A n symmetric with the same spectrum (shifted by a constant) as H An . ■ 

The long run distribution of the Markov chain described in Theorem 6 is given by the 
vector 7r = (tti, 7r 2 , ...), where 7Tj represents the "long time" probability of the state i G S, 
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satisfying 

(n+l) L (n+l) L 

E (QA n )af37T a = 0, Vf3eS, J2 ^ = 1- (73) 
a=l a=l 

However as this Markov chain is reducible, the solution of the equation ([F3|) is not unique 
but depends on the initial conditions. From (|35|) we see that 

(ti') a/3 = (h")p a = 0, V/?, a = l(n + l) + l + l, I = 0,1, ...,n. 

Hence from ( [72] ) we get 

(Q An) a/3 = (QA n )f3a = 0, a G S , V/3, 

with So as in (|68|). Therefore if this Markov chain is initially in a given state a G So, it 
will remain in that state. 

The states (3 $ So form a closed subset of S. From (]35f) and (|72|) one also learns 
that the absolute value of all the nonzero elements of any column of the intensity matrix 
Qa u are equal. Let S' be a closed subset of S with / elements. If the Markov chain is 
initially in the closed set S', then it will remain in 5" and the long run distribution is 
7r = (tci,7T2, VT( n+1 )L+i), where iii = l/l for i 6 S" and 7Tj = if i ^ S'. 



5.4 Discrete and Continuous-time Markov Chains Related to 
SU(2) Symmetric Integrable Ladder Models 

To every site on the z-th rung and #-th leg, 9 = 1,2, of the ladder we associate states 
described by the variable r it Q taking values and 1. The state space of this algebraic 
ladder is then finite and has a total of m — 2 2L states. 

[Theorem 7]. The following matrix 

1 L - x 
PsU{2) = 4(L-l)(18 + 4a + 46 + c) § (?4) 
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defines a stationary discrete-time S77(2)-symmetric integrable Markov ladder for a + 26 > 



16. The operator 




is g 


;iven 


by 


























( a i 


ft 2 


0,2 


a,-? 


ft -3 

a 


ft.i 


fti 


fli 


ftq 
o 


Q 4 


a,,i 


fli 


Clq 


fti 


fl„i 


a 4 ^ 






a 2 




ft IS 


ftR 


a 7 


aq 


ft* 


ft* 

u o 


ft* 


fto 

n y 


ft.i 


fli 


ft* 


fto 


fl„i 


a 4 






a 2 


ftfi 





ftfi 


ft* 


a 4 


fto 


fli 


fl,7 


ft* 

u O 


aq 


ft* 


ft* 


ft.i 


fto 


a 4 






a 2 


a 6 


a 6 




ft 8 


a 4 


fti 


ft 9 


ft 8 


«i 


a 4 


a, 9 


Cl- 


ft 8 


ft 8 


«3 






0-3 


ft 7 


a,* 


ft* 

o 


ftr. 


a 2 


ftfi 


ftfi 


fto 
y 


ft* 


ft.i 


a 4 


Cln 


ft* 


ft,! 


a 4 






a 4 


ftq 

o 


Q 4 


ft 1 


fl 2 


a l 




a 2 


a 4 


Clq 


fli 


a 4 


fti 


fti 
o 


fl„l 


a 4 






a 4 


ft* 


fto 

y 


ft 1 


ftfi 


ft 2 


ft^. 


ftfi 


fts 

o 




Clq 


a* 

o 


fti 


ft* 


fto 

y 


a 4 






a 4 


ft* 


ft 1 


fto 


ftfi 
v ' u 


ft 2 


ftfi 


11 


ft.i 


ft* 

v ' o 


ft.i 


fto 


ft* 


ft 7 


fl* 

u o 


a 3 




n" = 


































(75) 




^3 


ft* 

o 


ft 7 




fto 

.1 


a 4 


ft* 

o 


fli 


dec 



ftfi 




a 2 


ftfi 


flo 

y 


ft.i 


ft* 

o 


a 4 






a 4 


ftq 

v.* y 


ft* 


fli 


ft* 


ft 5 


ft 7 


ft* 


ftfi 





a 2 


ftfi 


ft.i 


fto 
y 


ft* 


a 4 






a 4 


a i 


Clq 


fli 


ft 4 


a 4 


ft 5 


C/ j 


ft? 


a 2 


a i 


0,2 


ft.i 


fli 


fti 


a 4 






a 4 


ft i 


u o 


fto 

y 


ft 4 


a 4 


I* ft 


I* y 


ftfi 


Or 




fl-2 


CIk 


ft* 


ft* 

o 


ft 7 


a 3 






a 3 


a 8 


«8 


(7,7 


u/y 


a 4 


fli 


«8 


flf) 


fti 


ft.i 


fl* 
u 8 


fl r 


Or 


Or 


ft2 






a 4 


a 9 


a 4 


a 8 


ft 8 


«3 


a 8 


ft 7 


a 4 


ft 9 


a.i 


a 8 


ftfi 


«5 


ft 6 


a 2 






a 4 


a 4 


a 9 


a 8 


a 4 


a 4 


a 9 


a 8 


fl 8 


ft 8 


«3 




ftfi 


a 6 


ft 5 


a 2 






\ «4 


fli 


a 4 




a 4 


a 4 


fli 


«3 


ai 


fli 


fli 


«3 


«2 


a 2 


ft 2 


ai j 




where a-i = 


66+a+46+4c, a 2 
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[Proof]. For the integrable ladder model (|^) with Hamiltonian H' — ^ii+i; the 

system remains integrable under (|60|) . It is straightforward to prove that 7i" = BTi'B -1 , 
where 



/-ll \ 

1 1/2 -1/2 1 

-1/2 -3/2 

y 1 -1 J 
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L-1 

Therefore the Hamiltonian systems H 1 and H" = ^ "H'li+i satisfy the relation (|60|) with 

i=i 

Bi = B, i = 1,2, ...,L. Hence H" is also S77(2)-symmetric and integrable with the same 
spectrum as H'. 

For a + 2b > 0, as the entries of Tt" are positive, H£p > 0, a,/3 = 1, 2, 2 2L . From 
(JT5P we also have Ea=iK/3 = 4(18 + 4a + 4o + c), V/? = 1,2, 16. By the definition (§7J) 
Psu{2) is the transition matrix of a stationary discrete-time S'?7(2)-symmetric integrable 
Markov processes. ■ 

The state space of this Markov processes associated with the stochastic matrix Psuti) 
is S — (1, 2, 2 2i ). Generally there is no closed subset C of the state space S such 
that {Psu(2))ij = for all i 6 C and j ^ C . However in certain parameter regions for 
a, b, c, from (|75|) one can see that there can exist such closed subsets C of S (the Markov 
processes are by definition reducible in these cases). From ([[^) it can also be seen that 
there exists no absorbing state for this Markov process. 

By using the results used in the proof of theorem 7, we have also the following inte- 
grable stationary continuous-time Markov process: 

[Theorem 8]. The matrix 

L-1 

Qsu(2) = H" - 4(L - 1)(18 + Aa + Ab + c) = - 4(18 + 4a + 46 + c)) M+1 (76) 

i=i 

is the intensity matrix of a stationary continuous-time Markov process. 

6 Conclusion and Remark 

Using the Casimir operators and coproduct operations of algebras, we have given a gen- 
eral way to construct square lattice models with a certain Lie or quantum Lie algebraic 
symmetry. As applications we discussed integrable A n symmetric chain models and SU(2) 
invariant ladder models. We have shown that the stochastic processes correspond to both 
A n symmetric integrable chain models and SU(2) invariant ladder models are exactly 
solvable stationary discrete-time (resp. continuous-time) Markov chains with the transi- 
tion matrices (resp. intensity matrices) which coincide with those of the corresponding 
integrable models. Other symmetric integrable lattice models ( e.g. with B n , C n , D n 
symmetry) and integrable Markov models can be investigated in a similar way. 
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